Electro-spinon in one-dimensional Mott insulator 
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The low-energy dynamical optical response of dimerized and undimerized spin liquid states in a 
one-dimensional charge transfer Mott insulator is theoretically studied. An exact analysis is given 
for the low-energy asymptotic behavior using conformal field theory for the undimerized state. In 
the dimerized state, the infrared absorption due to the bound state of two solitons, i.e, the breather 
mode, is predicted with an accurate estimate for its oscillator strength, offering a way to detect 
experimentally the excited singlet state. Effects of external magnetic fields are also discussed. 
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Recently, the optical activity of insulating magnets has 
been attracting intensive interests [H, 0, U 0] • Usually, 
the Mott insulator is considered to be inert in the lower 
energy region compared with the Mott gap. However, 
there are rich structures of the dielectric response of the 
insulating magnets in the infrared region as revealed by 
recent experiments in multiferroics [2, El, 0], which are 
explained partly by theories [f], . The central concept 
there is the electro- magnon, i.e., the combined wave of 
spin and polarization, which can be detected both by 
neutron scattering experiment and optical absorption. 
There are basically two mechanisms for the optical ac- 
tivity due to the electro- magnons. One is due to the 
spin-orbit interaction which is believed to be the 
origin of the multiferroic behavior in i?MnC>3 (i?=Tb, 
Dy) 0- The electro-magnon absorption has been first 
observed in these materials [2j. Although the spin cur- 
rent model based on the spin-orbit interaction Q explains 
well the static polarization in i?MnC>3, it turned out 
that this mechanism gives too weak oscillator strength 
for the infrared absorption. The other mechanism is 
the exchange striction where the polarization is given by 
P = H»j ft ' Sj m the absence of inversion symmetry 
at the center of the bond connecting the sites i and j [9( . 
This mechanism does not require the spin-orbit interac- 
tion, but usually gives rise to the two-magnon absorption. 
However, the recent finding is that it leads to the singlc- 
magnon absorption in non-collinear spin structures [6|. 
This can be easily seen by writing ft = (Si) + 5 Si, and 
consider the contribution to P linear in the fluctuation 
SSi from the ordered moment (Si). The linear terms can- 
cels when (Si) and (Sj) are collinear while they survive 
for the noncollinear structure. 

In this paper, we study theoretically the optical con- 
ductivity (j(u>) of the spin liquid state where the spin 
order is lost and the noncollinear spin structure is dy- 
namically generated. For that purpose, we propose that 
the (quasi) one-dimensional (ID) charge transfer Mott 
insulator is an ideal arena where optical spectra pro- 
vide a unique information on the quantum dynamics of 



the spins. In ID spin liquids, spinon (spin-i object) is 
the most fundamental particle [1 01 ] and the electro-spinon 
governs the infrared activity. Especially, in the undimer- 
ized state, the low-energy asymptotic behavior of u(ui) 
can be analyzed by conformal field theory (CFT). Fur- 
thermore, in the dimerized state, it is shown that the 
magnetically silent singlet excitations can be detected by 
optical means with a precise prediction of the oscillator 
strength. 

In the charge transfer compounds, the alternating 
stack of donor (D) and acceptor (A) molecules forms a 
chain. The simplified model for this system is the half- 
filled ionic Hubbard model 
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where A is the difference between the energy levels of D 
and A, U = to + g'ui is the transfer integral modulated 
by the molecular displacement Ui through the coupling 
constant g' , and <j>i is the electric potential at the molec- 
ular site i. The other notations are standard. Neglecting 
ti, the ground state is the ionic Mott insulator with one 
electron on each site when U > A. On the other hand, 
the ground state is neutral for U < A, i.e., all the odd 
sites are fully occupied while the even ones are empty. 

We are interested in the ionic state, and its effective 
spin Hamiltonian up to the linear order in u and (f> is 
given as follows: 

H eS = E J + 9Ui + vi-lTi^+i - 4>i)]Si ■ ft+i (2) 
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with J = AtlU/(U 2 - A 2 ), g = 8g't U/(U 2 - A 2 ), and 
7} = 8etQ[/A/(C/ 2 — A 2 ) 2 . Because of the lack of inversion 
symmetry at the center of the bond connecting i and 
we have the polarization along the chain, which is given 
by the derivative of H c g with respect to the electric field: 
dH cS 
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FIG. 1: Optical conductivity cr(a>) of ID Mott insulators. We take the unit J = C = a = 1. (a) cr(a>) of the XY spin chain 
for guo = 0.31 (dimerized) and that for guo = (undimerized; K = 2) at T = j-. Note that 77 = irr)i/v2 when K — 2. 
(b) Low-energy behavior of ct(uj) of the undimerized XXZ chain for K = 1, 3, and 2 at T = which is obtained from Eq. 
(£9j) by putting q — 0. The spinon velocity is given by u = n.Ja/2 — n/2 when K = 1. Due to the interaction effect, o"(w) is 
enhanced at low energies as K is decreased. Note that CFT description is valid only if u, T -C J = 1. In the limit of w 3> T, 
<t(oj) for if = 1 approaches to the constant value 2-Kr}\. (c) Low-energy behavior of (j(uj) of the dimerized Heisenberg chain 
with M = 0.2 and v = tt/2 at T = 0. Delta function peak is located at lo = \> r 3Mv. The s-s continuum starts at lo = 2Mv. 
Comparing (b) and (c), the continuum intensity is transferred to the breather peak. 



where a is the lattice spa cing 
as TTF-CA and TTF-BA " 



In real materials such 
12t the coupling to the 



molecular displacement Ui leads to the dimerization m = 
uq{— I) 1 accompanied by the singlet dimer state, i.e., the 
spin-Peierls instability occurs below a transition temper- 
ature Tsp. In this case, the uniform ferroelectric moment 
appears according to Eq. ([3]) . 

Let us now study the optical spectra in this system. 
From the Kubo formula, the optical conductivity ct(lo) is 
given by 



cr(w) 



C_ 

La 



wRe 



dte^- 5)t ([P(t),P(0)}}, (4) 



where L is the total number of sites in the chain and 
C = a/(hV u . c .) with the volume of the unit cell V n . c ,. 
Using the Jordan- Wigner transformation, the effective 
spin model in the absence of external fields is written as 
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off 
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where Jf = J a (l+(— l) l gua/J) for a =_L or z with J 1 - = 
J z = J. Here fi and // denote the spinless fermions, and 
n i = fili- The XY model with J z = is equivalent to 
the free-fermion model. It can be analyzed exactly and 
gives a useful starting point to analyze the realistic case 
J z = J^-. Supposing that the polarization operator for 
,V - is given by P = ^ (SfSf +1 + S?S* +1 ), 
the optical conductivity a(uj) is evaluated as 



a{u) = V^^-ujD^I-^) [f(-uj/2) - /( W /2)], (6) 

where oj\ — 2guQ is the spin gap induced by the dimeriza- 
tion, while L02 = 2 J is the bandwidth of the spinon dis- 
persion. The coefficient V is given by V = r} 2 Ca/[J 2 (l — 
lu 2 /lo 2 ) 2 ] and the Fermi distribution function is denoted 



by f(E) = (e E l T + with k B = H = 1. The behav- 
iors of a(ui) are shown in Fig. 1 (a). Note that in the 
undimerized state, <j(lo) oc u> 2 /T for lo <C T. 

This result for the non-interacting fermions should be 
modified by the interaction term, i.e., J z terms in Eq. 
([5]). Physically, this term introduces the attractive in- 
teraction between the particle and hole, which leads to 
the excitonic effect and enhances <j{ui) at lower energies. 
In the weak-coupling limit, this many-body effect can be 
studied using field theoretical methods as described be- 
low and several theoretical predictions will be made. 

It is well known that the low-energy physics of 
the undimerized Heisenberg chain is described by a 
Tomonaga-Luttinger (TL) model with a boson field & (x), 
which is equivalent to c = 1 CFT In addition, the 
polarization operator, i.e., the staggered dimerization op- 
erator, can be bosonized as 
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where e(x) = (— • S'i+i. Therefore, the continuum 
model for the dimerized Heisenberg chain Eq. @ is given 
by the sine-Gordon Hamiltonian [14j |: 



Hsg = ^ / 
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where H(x) is the canonical conjugate to &(x), which 
satisfies n(y)] = i8{x — y), the coupling constant 

fi oc guo, and K is the TL parameter. Here we have set 
a = 1 and the spinon velocity v = 1 . In the present nota- 
tion, the SU(2)-symmetric point (J 1 - = J z ) corresponds 
to A" = 1, while K > 1 to the easy- plane anisotropic 
case (J x > I J z \). In the high temperature region above 
Tsp, the lattice displacement uo is zero and hence we 
can set [i — in Eq. ([5]). Then the model is reduced 
to the c = 1 CFT. Since the conformal weight of e(x) is 
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(A, A) = (K/4:,K/4), the two-point correlation function 
behaves as (e(z)e(O)) ~ l/|z| with z — r + ix where 
t = it is the imaginary time. Using the finite-size scaling 
argument, its retarded Green's function at finite temper- 
ature T is obtained as 



G R (g,c) = -sin(^)(2^T)^ 2 FK(g, 
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whereM«, W ) = S(f-i^ff,l-fMf-i^,l-f) 
with the Euler beta function B(x, y) [15J. From Eqs. (j4|) 
and ((9]), we obtain a(ui) — —rj 2 CujImG R (0,uj) for several 
values of K as shown in Fig. 1 (b). The optical conduc- 
tivity is actually enhanced at low energies with a decrease 
in K, i.e., the increase in J z . Note that CFT description 
is valid only if uj and T are much smaller than v/a ~ J. 
For uj <C T, we find that a(uj) oc uj 2 T k ~ 3 . It is consistent 
with our earlier analysis of the XY model corresponding 
to K = 2. On the other hand, we find <j(u>) oc uj k ^ 1 for 
uj ^> T. It is remarkable that a(uj) is asymptotically con- 
stant at the Heisenberg point (K = 1). This constant is 
given by 2itri 2 Ca. It is possible to estimate the unknown 
coefficient rj\ from the experimentally observed a(oj). 

In the low temperature region below Tsp, /i in Eq. 
([5]) is nonzero and has a significant effect on the system. 
Fortunately, the model is still exactly solvable and one 
can apply powerful techniques such as the Bethe ansatz 
and form factor expansion [l6|. Let us concentrate on 
the realistic K = 1 case at T = 0, where the spectrum of 
Hsg consists of four quasi-particles, i.e., kink (s), anti- 
kink (s), and light (61) and heavy (62) breathers. The 
breathers are the bound states of the kinks and anti- 
kinks. Due to the hidden SU(2) symmetry, s, s, and 61 
form a spin-triplet with the mass gap M while 62 is a 
spin-singlet with the mass \[%M. The formula for the 
dimensionless gap M can be found in 
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where A ~ 3.041 and T(x) denotes the Euler gamma 
function. Physically, it corresponds to the spin-Peierls 
gap divided by hv/a = irJ/2. Suppose that the relation 
between the polarization density and the bosonic field is 
given by p(x) — rji cos(^/2n$(x)). Then the spontaneous 
polarization can be expressed by M as 



(p(x)) = (A/3)^ 2 Vl VM 



(11) 



using the vacuum expectation value of the vertex opera- 
tor e I/3 * 17j. Therefore, from the polarization measured 
in experiments, one can determine r\\\fM exactly. 

Next, we introduce the form factors and discuss the 
optical conductivity at T = 0. The Fock space of the 
sine-Gordon model can be constructed from the creation 
operators for the quasi-particles: \9i, fe"-^n)e 1 ,e 2 ...e n = 
Al(di)At 2 (d 2 )...Al n (9 n )\0). Here A\ k (e k ), the'ereation 



depends on the rapidity variable 9k which is related 
to the energy and momentum through the relations: 
Ek = M tk cosh0k and = M ek sinh^. The opera- 
tors A\.(9i) and A ej (9j) are neither bosons nor fcrmions 
and satisfy nontrivial commutation relations (lp| . In this 
Fock space, the form factor of cos(^/27^<^ , ) is defined as 



F cos (^...^) £l . 



cos(V27r$(O,O))|0i...0„> ei . 



The optical conductivity can be expressed in terms of the 
form factors as 



(12) 



To study the above expansion, it is useful to note that 
Hsg is invariant under the charge conjugation, C$C _1 = 
— <E>. From this symmetry and CA^C^ 1 = A\ 2 , the first 
nontrivial contribution stems from the heavy-breather 
form factor 14, 18, lij 2(J. The absolute square of this 
form factor is independent of the rapidity 9 and is given 
by 
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From Eqs. (fT2|) and (fT3|) . the optical conductivity due to 
the heavy breather is exactly given by 



I JT'cos 12 

a b2 (uj) = tttjICoj 1 b2 1 5(uj - VSM). 
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The next-leading contribution comes from the s-s con- 
tinuum (see Fig. 1 (c)). The threshold of this continuum 
is at oj — 2M. Therefore, the integrated intensity up to 
this threshold is given by 



Io M duja b2 {uj) 



■K7]lCv\F§° s \ 2 cx ri'fM 
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with the insertion of v and a. Combining this oscillator 
strength with the peak position located at ui = \/3M, one 
can exactly determine the coefficient r\\ experimentally 
At the same time, from the spontaneous polarization Eq. 
pip , one can also determine r\\ and compare it with 
that obtained from the oscillator strength. Therefore, 
the validity of the sine-Gordon theory can be tested by 
means of optical conductivity and spontaneous polariza- 
tion measurements. Let us now estimate the actual value 
of the integrated intensity for TTF-BA. This compound 
shows the spin-Peierls transition below Tsp = 53 K and 
the spontaneous polarization along 6-axis is measured as 
P s - 0.15/iC/cm 2 21]. Combining 
and K.c. ~ 8.4 x 4.3 x 11.6 A 3 \M 



with J ~ 170K 22] 



operator for = s, s, 61, or 62 with the mass M £fc , cillator strength as J^ M dojab2{u) ~ 10 13 Q 'cm x s 



we estimate the os- 
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FIG. 2: Optical conductivity <j(uj) of the field-induced TL 
liquid with magnetization M z = 0.1 for K — 1, |, and 2 
at T = We use the unit « = a = 1. The location of 
w c = 2irMzv/a is indicated. The imaginary part of G n (q,ui) 
is nonzero in the shaded region of the right figure. Note again 
that CFT description is valid only if u, T <C v = 1. 



which can be measured by infrared absorption spec- 
troscopy. 

Next, we shall consider the effect of the Zeeman cou- 
pling term 7? mag = —{JUbH J2i $f • As long as the mag- 
netic field h = gfisH is smaller than the spin-Peierls gap 
M, the singlet excitation (62) is insensitive to the field 
while the energies of the triplet excitations (s, s, and 61) 
are split. On the other hand, if h > M, the one species 
of the triplet particles becomes gapless and a uniform 
magnetization M z = (S?) appears. Its effective Hamil- 
tonian is equal to a TL model with a new boson field 
= — \/2ttM z x. The polarization density can 
be rewritten as p{x) = r\\ cos(v / 27r$(x) — 2ttM z x). Note 
that this TL-liquid state is adiabatically connected to 
that of the undimerized Heisenberg chain (uq = 0) under 
the magnetic field. The magnetization dependence of the 
TL parameter K in dimerized and undimerized Heisen- 
berg chains can be found in [24| and (25|, respectively. 

Applying CFT, we can again evaluate the correlation 
function of the polarization as in Eq. ([9]). The two-point 
correlator behaves as (p(z)p(O)) ~ cos(2irM z x)/\z\ K , and 
the retarded Green's function around the wavenumbcr 
q = is represented as G R (q,u) = G R (q — 2itM z ,uj) + 
G R (q + 2ttM z ,uj). From this, we can obtain <j(ui) = 
— rj\C uilmG R (0, ui) as shown in Fig. [2 Remarkably, 
in sufficiently low temperatures T <C v/a, u{lj) starts 
from a finite frequency around u c = 2irM z v/a. This is 
in sharp contrast to the gapless behavior in the zero- 
magnetization case [see Fig. ffl(b)]. 

Lastly, we compare the above results with those of Ra- 
man scattering j26| . The Raman spectrum is propor- 
tional to the dynamical structure factor of R cx JY Si ■ 
Si+i for the spin chain Eq. @. For the dimerized case, 
the Hamiltonian does not commute with R and hence its 
spectrum is equivalent to the dynamical structure factor 
of R cx J2i(~ l)'^ ' Si+i. Therefore, the Raman scat- 
tering can also provide a spectrum with a singlet peak 
as has indeed been observed in CuGeC>3 



the intensity of the Raman spectrum, in general, while 
we can evaluate that of the absorption with consider- 
able accuracy as in Eq. (|14| . Furthermore, the infrared 
absorption can detect gapless spinon excitations in the 
undimerized case where any contribution from spinons 
is absent in the Raman scattering since the operator R 
commutes with H c g with m = 0. 

In summary, we have studied the optical activity of the 
electro-spinons in the undimerized and dimerized quan- 
tum spin chains with S = \ derived from the ionic Hub- 
bard model. The scaling form of the low-energy <j(u>) is 
predicted by applying the conformal field theory in the 
undimerized case, and the exact estimate of the oscilla- 
tor strength due to the singlet breather mode is given. 
These theoretical predictions can be tested in organic 
charge transfer Mott insulators. 
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we should note that it is difficult to theoretically predict 
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